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Abstract 

It is shown that the known method to solve the Dirac equation by means of the squaring method, 
when relying on the scalar function of the form $ = e _1£t e* fcix e lfc2!/ sin(/c2+a) leads to a 4-dimensional 
space of the Dirac solutions. It is shown that so constructed basis is equivalent to the space of the Dirac 
states relied on the use of quantum number fci,fc2,±fc and helicity operator; linear transformations 
relating these two spaces are found. Application of the squaring method substantially depends on the 
choice of representation for the Dirac matrices, some features of this are considered. Peculiarities of 
applying the squaring method in Majorana representation are investigated as well. The constructed 
bases are relevant to describe the Casimir effect for Dirac and Weyl fields in the domain restricted 
by two planes. 
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1 Introduction 

In connection with the Casimir effect [1[ for the Dirac field in the domain restricted by two planes, 
of special interest are solutions of the Dirac equation, which have vanishing the third projection of 
the conserved current J z on the plane boundaries. Such solutions are reachable when considering 4- 
dimensional space 


{\k} = {'&kiM,ka=k,<r ® ^fci,fc 2 ,fc 3 =—fe,cr}, (1-1) 

that is the basis of four solutions with opposite signs of the third projection of momentum +^3 and —^ 3 ; 
<7 is referred to polarization of the states. 

It is shown that the known method to solve the Dirac equation trough squaring method [7] (elaboration 
of such a method to the case of electromagnetic field see in 001 ), when relying on the scalar function of 
the following form 

<f> = e - iet e iklX e ik2V sm{kz + a) =► {$} = {$ 1 , tf 2) $ 3 , $ 4 } ( 1 . 2 ) 

leads to a 4-dimensional space of the Dirac solutions. It is shown that so constructed basis ED is 
equivalent to the space of the Dirac states (El); linear transformations relating these two spaces are 
found. Different values of parameter a in ED determine only different bases in the same linear space 

ED- 

Application of the squaring method substantially depends on the choice of representation for the Dirac 
matrices, some features of this circumstance are considered. 
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Peculiarities of applying the squaring method in Majorana representation are investigated. It is 
shown that constructed on the base of scalar functions cos (et — kx) and — i sin(et— kx), two 4-dimensional 
sets of real and imaginary solutions of the Majorana wave equation cannot be related by any linear 
transformation. 

General conditions for vanishing third projection of the current J z at the boundaries of the domain 
between two parallel planes are formulated: firstly, on the basis of plane spinor waves ED; and secondly, 
on the squared basis (11.21) . In both cases, these conditions reduce to a linear homogeneous algebraic 
systems which leads to a 4-th order algebraic equation, the roots of which are e 2lfea , where a is a half¬ 
distance between the planes, and k stands for the third projection of the Dirac particle momentum. Each 
solution of this equation, with represent a complex number of the unit length, provides us with a certain 
rule for quantization of the third projection of k. Explicit forms of these algebraic equations are different, 
however their roots must be the same. 

Conditions for vanishing the current J z for Weyl neutrino field on the boundaries of the domain 
between two planes are examined, the problem reduce to a 2-nd order algebraic equation. 

Covariantization of the conditions for vanishing J z is performed. 


2 Squaring method 

Let us start with a special solution of the Klein-Fock-Gordon equation with one formal change in 
4> = e ~ Kt + kx : namely, we will change the factor e +lksZ by the real-valued factor 

sin(kz + 7 ) = sin<^ ; ( 2 . 1 ) 

for brevity we use notation k = k^. Thus, we start with 

(■ H a da - M)(i r y a d a + M) = (- d 2 t + d,dj - M 2 ) , 

$ = e -itt e ik 1Xe ik 2 y sin ^ ) e 2 _ fc 2 _ kj _ k 2 - M 2 = 0 . (2.2) 

The 4x4- matrix of the four columns-solutions of the Dirac equation is constructed in accordance with 
the following rule 

{'Ll, ^ 2 , 'Ls, ^ 4 } = (h a d t + h j dj + M)$ . (2.3) 


In spinor basis, at the choice (12.21) . the matrix (12.31) will be 


vj> 1 _ e ~i£t £ ikix e ik 2 y 

M siny: 

0 

(esin<^ + ik cosip) 

— (fci + 2 ^ 2 ) sin ip 

4/2 = e~ iet e iklX e ik2V 

0 

M sinyj 

— (fci — ik 2 ) sin ip 
(esintp — ik cos ip) 

(]> 3 _ e ~itt e ik!X e ik 2 y 

(esin<p — ik cos ip) 
(fci + ik2) sin^ 

M siny; 

0 

, 4/ 4 = e ~ iet e ik ^ x e^v 

(fci — ik 2 ) sin p 
(esint^ + ik cosip) 

0 

M sintp 


To decide on the linear dependence of 'L j (or not), one should examine the following relation 

ai’Li + q > 2 41 2 + U3'I f 3 + < 244/4 = 0 , 

so we get the linear homogeneous system with respect to ai, <22,(13, <24: 

a\M siniy? + <2,3(6 sin p — ik cosip) + 04(^1 — 2^2) sintp = 0 , 

<22 M sin<^ + 03(^1 + 1^2) sine/? + <24(6 sin + ik cosip) = 0 , 
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ai(esint ( 5 + ikcostp ) — 02 (fci — *^ 2 ) siny> + 03 M sin 73 = 0 , 

—a\{ki + ik 2 ) sin </3 + a 2 {ts\ntp — ikcostp) + a 4 M sin</j = 0 . (2-5) 

Its determinant turns to be det{'I'} = k 4 . This means that four solutions (12.41) of the Dirac equation (at 
any a in the tp = kz + a) are linearly independent: they determine a 4-dimensional space. 

Let us specify two possibilities for 7 in (12.11) : 


7 = 0 , tp = kz, sin tp = sin kz, cos tp = cos kz , 


V ]> 1 _ e ~ iet ,gikix e ik 2 y 

M sin kz 

0 

(e sin kz + ik cos kz) 

— (fci + ik 2 ) sin kz 

vJ / 2 _ g-iet e ik x x e ik 2 y 

0 

M sin kz 

— (fci — ik 2 ) sin kz 
(e sin kz — ik cos kz) 

_ g-ietgifcix e ik 2 y 

(e sin kz — ik cos kz) 
(fci + ik 2 ) sin kz 

M sin kz 

0 

, T 4 = e - i£t e ifclX e ifc2y 

(fci — ik 2 ) sin kz 
(e sin kz + ik cos kz) 

0 

M sin kz 


7= 2 ’ 

V]/^ _ e ~iet e ikix e ik 2 y 


tp' = kz — 


2’ 

—M cos kz 

0 

(—e cos kz + ik sin kz) 
(fci + ik 2 ) cos kz 


sin tp' = — cos kz, cos tp' = sin fez , 


$' (2) = e -izt e ikix e ik 2 y 


0 


—M cos kz 
(k\ — ik 2 ) cos kz 
(—e cos kz — ik sin kz) 



(—e cos kz — ik sin kz) 


— (fci — ik 2 ) cos kz 

^ = g-iet gifei x e ik 2 y 

— (fci + ik 2 ) cos kz 
—M cos kz 

_ g— ietgjfciXgife 2 y 

(—e cos kz + ik sin kz) 

0 


0 


—M cos kz 


We readily find expressions for the following combinations (the total factor e let e lk i x e lk 2 y i s omitted): 


- T( + 'bl/ , = e ikz 


M 


0 

0 

e — fc 

, -%+i^ 2 = e ikz 

M 

-(fci - ik 2 ) 

— (fci + ik 2 ) 


e + k 



e + fc 


(fci - ik 2 ) 

^ + i4 - 3 = e ikz 

(fci + ik 2 ) 
M 

, -% + M 4 = e ikz 

e — fc 

0 


0 


M 


( 2 . 8 ) 


Note that solutions (12.81) exactly coincide with those obtained by applying the squaring method when 
one starts with a scalar function as the ordinary plane wave: 


U = 


^-iet-\-ikix-\-ik 2 y-\-ikz 


M 

0 

(e-k) 

— (fci + ik 2 ) 


0 (e + k) 

M (fci + ik 2 ) 

— (k\ — ik 2 ) M 

(e + fc) 0 


(fci - ik 2 ) 
(e - k) 

0 

M 


(2.9) 


The rank of the matrix in (12.91) equals to 2. Therefore, among four solutions of the Dirac equations 
given by (EH) only two of them are linearly independent. For definiteness, let us chose solutions and 
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U( 2 ) : 


Ui = e 



M 


0 

iet-Hfci ai-Hfo y-\-ikz 

0 

tt _ — £et-t-ifci£C+ifc22/+&fc2: 

M 


e — k 

, U2 — " 

-(ki - ik 2 ) 


-{k\ + ik 2 ) 


e + k 


It is readily checked that remaining solutions U 3 , 1/4 are expressed through U\, U 2 as follows 
U 3 = -jj [ (e + k) Ui + (fci + ik 2 ) U 2 ] , £4 = [ (fci - ifc 2 ) Ui + (e - fc) C/ 2 ] . 

Also, instead of (I2.8[> . one can construct other combinations 



M 


0 

T) - iTi = e“ ifcz 

0 

e + k 

— (fci + ik 2 ) 

, - *$ 2 = 

M 

~{ki - ik 2 ) 
e — k 


% -W 3 = e 

e — k 
(ki + ifc 2 ) 
M 

, -^4 - i ^ 4 = e ife2 

(fci - ifc 2 ) 
e + k 

0 


0 


M 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


These coincide with those obtained by applying the squaring method to a scalar function as the ordinary 


(2.13) 


plane wave with the only change k 

into —k: 





M 

0 

(e- fc) 

(fci — */c 2 

jjt _ ^-iet-\-ik\x-\-ik' 2 y—ikz 

0 

M 

(fci + ife 2 ) 

(e + fc) 


(e + k) 

-(ki - ik 2 ) 

M 

0 


— (fci + ik 2 ) 

(e-fc) 

0 

M 


Again, the rank of the matrix in (12.1311 equals 2, so only two solutions are independent: 



M 


0 

jjf _ g— iet+2/cirr-|-i/c2 2/ — 

0 

e + fc 

— (fci + ik 2 ) 

tt/ -iet-\-ik\x-\-ik 2 y—ikz 

1 u (2) ~ e 

M 

~(ki - ik 2 ) 
e — k 


Solutions 1 /( 3 ) and Z7( 4 ) are constructed through £/(i),£/( 2 ) in accordance with the rules 


(2.14) 


u 'w ~ M L 


(e — k) U't + (fci + ik 2 ) C/( 2 ) 


= M L 


(fci - ik 2 ) U ^ + (e + fc) U[. 


( 2 ) 


(2.15) 


3. Solutions in the basis of momentum 4-vector and helicity 


Solutions of the Dirac equation in Cartesian coordinates can be searched in the form 

T £)fcljfe2 = e~ iet e iklX e ik2V e ik3Z 

With the use of spinor basis for Dirac matrices, one gets the linear system for /,: 


fi 

h 

h 

h 


e /3 + & 1/4 - ik 2 fi + k 3 f 3 - Mfi = 0 , 

e / 4 + kif 3 + ik 2 f 3 - k 3 f 4 - Mf 2 = 0 , 

e/i - k x f 2 + ik 2 f 2 - & 3/1 - M/a = 0 , 

e /2 - kifi - ik 2 fi + k 3 f 2 - M/ 4 = 0 . 


(2.16) 


(2.17) 
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Let us diagonalize the known helicity operator E = <JjPj. With the substitution (12.161) . from eigenvalue 
equation E $ = p $ we arrive at 


kih - *fc 2 / 2 + k 3 fi = pfi, 
kifi + ik 2 fi - k 3 f 2 = p /2 , 
fci/ 4 - ik 2 f 4 + k 3 f 3 = pf 3 , 

kif 3 + ik 2 f 3 - k 3 f 4 = pf 4 . (2.18) 

Considering eqs. (12.171) and (12.181) jointly, we obtain the system (note that p 2 = e 2 — M 2 ) 

e /3 + p /3 - Mf 4 = 0 , e/ 4 + p/4 - M f 2 = 0 , 
e/i - p/i - Mf 3 = 0 , ef 2 - pf 2 - M f 4 = 0 . 


This system results in two values for p and corresponding restrictions on /,: 


h 


e-p 

M 


/1, 


U 


e~P 

M 


h- 


Allowing for (12.19|1 . from (12.181) one gets the system for f\ and f 2 . 


(k 3 - p)/i + ( h - ik 2 ) f 2 = 0 , 
(k 3 + p)/ 2 - ( fci + ik 2 )/i = 0 . 


Further we obtain 


h 


i 

iki + k 2 


(k 3 - p)/i 


k 4 + ik 2 
k 3 +p 


fi- 


(2.19) 


( 2 . 20 ) 


Thus, two independent solutions (we will mark them by a and /3) at the fixed k are (further let f± = 1 
and p = +Ve 2 - M 2 ) 


(«)> 


e - P , , , h+ik 2 

^ J J 1 / 2 7 ^ ? 


M 




fci + 

/1 = 1, / 2 = r--=t; 

k 3 -p 


these may be presented in a shorter form 



1 


1 

V]> ( J _ e -iet e ik 1 x e ik 2 y e ik 3 z 

s 

a 

, T (/ 3 ) = e ~ iet e ik ix e ik 2 y e ik 3 z 

t 

(3 


as 


fit 


( 2 . 21 ) 


( 2 . 22 ) 


3 Relationships between two bases: those obtained from squar¬ 
ing method and the momentum-helicity solutions 

Four types of the different solutions of the Dirac equation can be constructed be the method of separation 
of variables (we take k 3 = k and k 3 = —k) 



1 


1 


ki+ik2 


ki+ik2 

= « < ia) (k)=e ik * 

k+p 

a 

k\+ik 2 

U k+p 

,$2 = 'S' ia )(-k)=e- ikz 

— k+p 

a 

k±+ik 2 
a -k+p 
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1 


1 


fci +ifc 2 


ki -\-1k2 

$ 3 = ^)(fc) = e ifez 

k—p 

P 

g k 1 +ik 2 
^ k—p 

,$4 = *( 0 )(-k) = e- ikz 

—k—p 

P 

g k\+ik 2 
^ —k—p 


and they must be related to four squaring solutions in Section 2. 

By physical reason, we should expect existence of the following linear expansions 


Ui = a$i + &<1>3, 

Evidently, if a, 6, c, d are known, one can derive 


£1 2 = c$i + d$ 3 . 


£1(3) — ^(1) + (fo + *^2) £1(2)] 

= —y [(e -t- fc)(a$i + 6$ 3 ) + (Aq + i/c 2 )(c$i + d$ 3 )] ; 


( 3 . 1 ) 


(3.2) 


(3.3) 


£ 1 ( 4 ) = — [(k 1 - ifc 2 ) C/ ( i) + (e- k)\U{ 2 )} 


1 


= -^[(&i _ *fc 2 )(a$i + 6$ 3 ) + (e - fc)(c$i + d$ 3 )] . 
Analogously, there must exist expansions 

JJ[ = a'$ 2 + 6'$ 4 , U 2 = c'$ 2 + 4 , 

with the help of which one can derive 

£ 1 ( 3 ) = ^[( e — k)(a'$ 2 + ^$ 4 ) + (Aq + *fc 2 ) d( 2 )] j 
£ 1 ( 4 ) = [(Aq — ifc 2 )(o , $2 + ^$ 4 ) + (e + fc)(c , $ 2 + d$ 4 )] ■ 
In turn, after that, one can derive the next expansions by the formulas 

Let us consider the first relation in (13.21) . U\ = a$i + £$3 ; explicitly it reads 


M 

0 

e — k 

— (Aq + *fc 2 ) 



1 


1 


ki+ik 2 


ki+ik 2 

= a 

k+p 

+ 6 

k—p 


a 

P 


ki~\~ik 2 

k+p 


g ki +ik 2 
^ k—p 


from whence it follows 

M = a + 6, 0 = 


k + p k — p 
From the first and second equations we get 


e — k = aa + 6/3, 


-1 = 


6/3 

.-b - - 

k + p k — p 


a=^(fc+p), 6=-^(fc-p); 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 


(3.9) 


(3.10) 


(3.11) 


two remaining equations turn to be identities. 
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Now, let us consider the second equation in (13.21) : U 2 = cd> 1 + 3 ; it reads explicitly 


0 


1 


1 

M 

-(h - ik 2 ) 


fci +ifc 2 


k-\ +ik2 

= c 

fc+P 

a 

+ d 

k—p 

P 

e + k 


fci+jfc 2 

U k+p 


0 k\+ik 2 
r k-p 


or 


0 = c+ d, 

— (fci — ik 2 ) = ca + d/3 , 


ki + ik 2 , M + ifc 2 
M = c— -1- d- 


k = 


k + p k — p 

ca(ki + ik 2 ) d/3(ki + ik 2 ) 


P 


k — p 


From the first and second equations, it follows 


d = —c , 



ik 2 ) ; 


(3.12) 


(3.13) 


two remaining ones are identities. Thus, we have obtained 


M M 

Ui = — [(p + fc)$ 1 + (p - fc)<h 3 ] , U 2 = — [(fci - ik 2 )<$> 1 - (fci - ik 2 )$ 3 ]; 


and further 


Mp+k p-k 

U 3 = -[- -$1 + $3] 


2p a 


P 


M k\ —ik 2 k 3 -ik 2 

Ui = —- $1 -——$ 3 ] 


2 p a 


P 


the identities e — p = M/(3, e + p = M/a should be remembered. 

Now, let us consider the equation U[ = a'<& 2 + b'$ 4 ; it reads explicitly as 


M 


1 


1 

0 


ki~\~ik2 


k-\ +ik2 

= a' 

-k+p 

+ b' 

—k—p 

e + k 

a 

P 

— (fci + ik 2 ) 


fcl+ifc 2 
a -k+p 


a k\+ik 2 
r -k-p 


(3.14) 


(3.15) 


(3.16) 


what differs formally from (13.91) only in the change fc —>■ — fc a presence of primes at variables, so the 
result can be written at once 


Consider the equation U 2 


a' = — (-k+p) , b' 
2 P 

c 7 <&2 + d'$i ; it reads 


M ( + ^ 


0 


1 


1 

M 

-(ki - ik 2 ) 


ki~\~ik2 


k\ +ik2 

= c 

—k+p 

a 

+ d' 

—k—p 

P 

e — k 


n ki+ik 2 

U -k+p 


q k^ +ik 2 

F -k-p 


which differs from (13.121) only in notation, so its solution looks d! = —c 7 , c! = ||(fci — ik 2 ). 
Thus, we have derived decompositions 


U[ = ^ [{p - fc )$ 2 + {p + fc)4> 4 ] , U' 2 = ^ [(h - ik 2 )$ 2 - (fci - *fc 2 )$ 4 ] • 


U' z = — 

3 2 p 


M (p — k 


$2 + 


P- 


-< 1>4 


t/ 7 -^ 


fci — ik 2 , fci — *fc 2 * 

-$2--- $ 4 


a 


P 


(3.17) 


(3.18) 
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and repeat (13.141) . (13.151) : 


M M 

U\ = —— [(p + fc)4>i + (p - , U 2 = — [(fcl - *fc 2 )$l - (fcl - *fc 2 )$3 

2p 2p 


M f p + k P-k 

u 3 = 7 r \ -+ —5— $ 3 

2 p \ a p 

Next, relying on f the formulas 

- U[), ^ 1 = ^U 1 + U[) 

^3 = 2 i ^ 3 ~ Uz)i ^3 = “2^ 3 + ^3) ) 
we derive decompositions 

4 h = ain^n , 

where the involved matrices are given by 


M (/cj - ik 2 h - ik 2 

C / 4 = — - $1 - - - $ 3 

2 p \ a p 


* 2 = ^(U 2 -U^, % = ~(U 2 + U[) , 

4-4 = 4(^4 -Ui), % = - l -{U i + U' A ) 


2 i 

4 '[ = a' ln <f> r , 



(P + k) 

-(p-k) 

(p-k) 

-(p + k) 

M 

(ki - ik 2 ) 

-(k\ - ik 2 ) 

-(ki - ik 2 ) 

(h - ik 2 ) 

4p 

a~ 1 (p + k) 

—a _1 (p — k) 

(3~\p-k) 

-j3~ 1 (p+ k) 


a~ 1 (k\ — ik 2 ) 

—a~ 1 (ki — ik 2 ) 

-/3 -1 (fci - ik 2 ) 

(3~ l (ki - ik 2 ) 


(p + k) 

(p - k) 

(p-k) 

(p + k) 

M 

(ki - ik 2 ) 

(ki - ik 2 ) 

-(ki - ik 2 ) 

-(ki - ik 2 ) 

4p 

a~ x (p + k) 

a _1 (p — k) 

p- x (p-k) 

(3 ~ l (p + k) 


a~ 1 (ki — ik 2 ) 

a~ 1 (k 4 — ik 2 ) 

-/3~ 1 (k 1 - ik 2 ) 

-p~ 1 (k 1 - ik 2 ) 


From those we can construct a new matrix 

M 


Sii — a-i-i + idij — 


2 P 


(p+k) 0 

(fci — ik 2 ) 0 
a -1 (p+k) 0 


(p — k) 0 
— (hi — ik 2 ) 0 

/3~ 1 (p—k) 0 


L (fci — ik 2 ) 0 —j3 1 (k\ — ik 2 ) 0 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(it corresponds to the use in the factor e lkz in the initial scalar substitution for 4>) The rank of this matrix 
equals 2 , and it is responsible for transformation 


u 1 : u 2 ,u 3 ,u 4 


4*1,4> 


3 • 


Analogously, we have another combination (it corresponds to the use in the factor e lkz in the initial 
scalar substitution for 4>) 


S ij - a ij - 2p 


0 (p — k) 0 (p + k) 

0 (ki — ik 2 ) 0 — (k\ — ik 2 ) 

0 a~ 1 (p—k) 0 ft -1 (p + k) 

0 a -1 (fci — ik 2 ) 0 —/3 _1 (fci — ik 2 ) 


(3.24) 


the rank of this matrix equals 2 , and it corresponds to the transformation 

s 


U (, U 2 . [/.(, 1/4 


4 > 2 , $4 . 


Because the determinants do not vanish 




















inverse matrices exist and have the form 


2 ip 


^ 1 kM(a - p) 


a 

a 

p 

p 


a.(k—p) 
ki—iko 
a(k+p) 
ki — iko 

P(k+p) 

k\—iko, 

P{k-p) 

k i —ik<2 


(k-p) 

ki —ik2 
(k+P ) 

ki —ik2 
( k+p ) 
k\ —ik2 
{k-p) 
k i —ik2 


(3.26) 


[(<*%] = 


2 ip 


kM(a — P) 


u(k-p) 
ki —ik2 
a (k+p) 
ki—iko 
_n (3 (k+p) 

^ k-\—ik2 

a P(k-P) 

k' k\ —ifco 


l 

-l 

l 

-l 


{k-p) 

ki —ik2 
(k+p) 

k\ —ik2 
(k+p) 

k\ —ik2 
(k-p) 

ki —ik2 


(3.27) 


Evidently, we have inverse transformations relating 'I > j and : 


a 'ik ( a )kl = 


2 k(a — P) 


[ajk+ i )fei] = 


k(a — P) 


% = a W l> k = 

[ a jk ( a l )kl 

* k , 



(p + k)a 
( k\ — ik 2 )a 

( k+p)P 

(p—k) 2 a 
ki —ik2 

— {k + p)a 
{p 2 ~k 2 )P 

~(p-k) 

(h - ik 2 ) 

p—k 

2 7 2 

p—k 

ki —ik2 

-(k+p) 

{p+k) 2 

, (3.28) 

(k 1 -ik 2 )P 

a 

(k\—ik2)oi. 

(p—k)/3 

a 

p 

k\ —ik2 

p 

{k\—ik 2 )P 

p—k 

~T 


^ - j — djk^k 

\^jk )/cZ 

K,, 



— (p + k)a 
(k\ — ik 2 )a 
~(k+p)P 

(p—k) 2 a. 
k 1 — ik2 

(k +p)a 
(p 2 ~k 2 )P 

(P - k) 

(ki - ik 2 ) 

p—k 

P 2 -k 2 

k\ — ik2 
(k+p) 
{p+k) 2 

(3.29) 

(k 1 -ik 2 )p 

(ki —ik2)ot 

{p—k)P 

p 

k 1 —ik2 

(ki—ik2)fi 
p—k 



The following conclusion can be given: the choice of an initial phase 7 in the function sin(fcz + 7 ) 
does not influence on the whole structure of the space of solutions - it only determines a specific basis in 
the same space: 


sin (kz + 7) = cos 7 [sin fez] — sin 7 [— cos kx\ , 
sin(fcz + 7 ) => d'j = cos 7 ’I'j — sin 7 tp' . 


(3.30) 


4 Dependence of squared solutions on the choice of Dirac ma¬ 
trices, standard basis 

Let us follows two representations for Dirac matrices 

7 a , r Q = s+s~ 1 . 

Having applied the squaring method in both bases, we will obtain two sets of solutions: 

(n°a a + M)$ = {^ 1 ,^ 2 ,T 3 ,^ 4 }, (ir a a a + M)$ = {^,^ 2 ,^3,^ 4 } ■ (4.i) 

These two sets, or two matrices, are linked to each other according to the formula 

S'i'f'i, ^2, ^ 3 ,1'd-S ' -1 = {<Pi, ■ (4.2) 
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This relation tells that if at a fixed initial scalar function $ we have a matrix of squared solutions 
with rank 2, the corresponding matrix of solutions in any other Dirac basis will have also the rank 2. 
Analogously, if we have the matrix of solutions of the rank 4 in one Dirac basis, we will have in any other 
Dirac basis a matrix of solutions of the same rank 4. 

For instance, let us specify squared solutions in the commonly used standard representation for Dirac 
matrices: 


0 


I 0 


0 cd 

7 


0 -I 

, i 1 = 

-u l 0 

Taking into account the relation 

{i^da + M) = 

id t + M 

0 

0 

idt + M 

id 3 

id 1 — d 2 


-id 3 

—idi - d 2 

- id t + M 


- 

id\ + d 2 

id 3 

0 


id i + d 2 
-td :i 
0 

-id t + M 


and choosing 


$ = e let e zklX e lk2V sin<£> , siniys = kz + 7 , 


we get an explicit form for the matrix of solutions 


m 


(e + M ) siniys 

0 

—ik cos p 
(fcq + ik 2 ) sin:/? 


0 

(e + M) sin^ 
(fci — ik 2 ) sin p 
+ik cosy) 


+ik cos p 

(—ki — ik 2 ) sin</? 

(—e + M) sin^ 

0 


(—hi + ik 2 ) sin p 
—ik cos p 

0 

(—e + M ) siny) 


Thus, we construct four different solutions of the Dirac equation: 


(4.3) 


(e + M) sin ip 


0 

0 

, w 2 = 

(e + M) sin<^ 

—ik cos ip 

(k\ — ik 2 ) sin^ 

{k\ + ik 2 ) sin <p 


+ik cos (p 



+ik cos p 


{—k\ + ik 2 ) sin p 

W 3 = 

{—k\ — ik 2 ) sin:/? 
(—e + M ) sin^j 

, w 4 = 

—ik cos p 

0 


0 


(—e + M ) sinyj 


Due to statement after (14.21) . they are linearly independent. 

Two choices of the phase 7 lead to two sets of solutions: 

p = kz, sin p = sin kz, cos p = cos kz , 


(e + M) sin 


0 

0 

, w 2 = 

(e + M) sinfcz 

—ik cos fcz 

( k\ — ik 2 ) sin kz 

(fci + ik 2 ) sin fcz 


+ik cos kz 



+ik cos kz 


(—ki + ik 2 ) sin kz 

W 3 = 

{—k\ — ik 2 ) sin kz 
(—e + M) sin fcz 

, w 4 = 

—ik cos kz 

0 


0 


(—e + M ) sinfcz 



sin tp' = — cos kz, cos p = sin kz , 


(4.4) 


(4.5) 
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(4.6) 



— (e + M) cos kz 


0 

w[ = 

0 

—ik sin kz 

, w^ = 

— (e + M) cos kz 
— (hi — ik 2 ) cos kz 


— {ki + ik 2 ) cos kz 


+ik sin kz 


+iksinkz 


(k\ — ik 2 ) cos kz 

(fci + ik 2 ) cos kz 

, wi = 

—ifcsin kz 

(e — M) cos kz 

0 

0 


(e — M) cos kz 


From these, one easily constructs linear combinations referring to application of the squaring method 
for two differen choices of a scalar function: 


e -iet e i kl x e ik 2 y [. _ W I + iWj ] = ( i7 « and Q a + e -itt e ik lXe ik 2 y e +ikz , 

e -iet e ikix e ik 2 y y_ w , _ = (^a + e ~ict £ ik lX£ ik 2 y £ -ikz _ 


(4.7) 


Evidently, one can perform the analysis like given after (12.91) (12.151) with some minor technical alter¬ 
ations. 


5 Squaring method and Majorana fermion 

Special interest for the method has any Majorana basis for Dirac matrices. Let us specify one of them 
as follows 


'I'm = A^ 


spinor ) 


A = ± dL , A - 1 = 1+72 


V 2 ’ 


V2 ’ 


r a M = A^A - 1 


explicitly these matrices read 


7m = +7°7 2 


7m =7 


2 


0 

—i 

0 

0 


—i 

0 

0 

0 

i 

0 

0 

0 

1 , 1 2 

0 

i 

0 

0 

0 

0 

0 

i 

, 7m = +7 7 = 

0 

0 

—i 

0 

0 

0 

—i 

0 


0 

0 

0 

i 


0 

0 

0 

i 


0 

i 

0 

0 

0 

0 

—i 

0 

3 , 3 2 

i 

0 

0 

0 

0 

—i 

0 

0 

, 7m = +7 7 = 

0 

0 

0 

i 

i 

0 

0 

0 


0 

0 

i 

0 


(5.1) 


First, let us construct real and pure imaginary solutions of the Dirac equation, staring from the Dirac 
solutions with n momentum-helicity quantum numbers. To this end, it is enough to translate the Dirac 
Diane waves (12.221) of the tvDes a. 3 from SDinor basis to Maiorana one ED, and after that we are to 
separate real and imaginary parts of these solutions. In this way we will construct wave functions for 
Majorana particles with different charge parities. 

Thus, we start with the plane waves of the form 


’J'(a) = e~ iet e ikS 


U, _ —ietikx 

V((3) — e e 


1 

s 

a 

as 

1 

t 

13 
131 


e — p _ ki + ik 2 
M S k3+p 


_ e+ p ki + ik 2 

M 1 &3 — p 


(5.2) 
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With the help of the matrix A ( 15 . 11 ) we translate them to the Majorana representation: 



(1 + ias) 


(1 + i/ 3 t) 

g—ikx 

(s — ia) 

g—ikx 

’ = V2 

(t - iP) 

72 

—i(s + ia) 
i( 1 — ias) 

—i(t + if}) 
i( 1 — ipt) 



(1 — ias*) 


a -m*) 

p -\-ikx 

<Tf * — 

(s* + ia) 
i(s* — ia) 

g-\-ikx 

{t* + W) 
i(t* — i/ 3 ) 


—i( 1 + ias*) 


-i(l + ipt*) 


( 5 . 3 ) 


( 5 . 4 ) 


where kx = et — k x. 


Then, from wave functions ( 15 . 3 |) . one separates real and imaginary parts: 


R M{a) = ^(^M(a) + ^M(a)) 


and 


r 

1 + ia(s — s*)/2 


a(s + s*)/ 2 


cos kx 

(s + s*)/2 
a — i(s — s*)/ 2 

+ sin kx 

—i(s — s*)/2 — a 



— (s + s*)/2 


{ 

a(s + s*)/2 


1 — ia(s — s*)/2 

> 


1 

72 


lM( a ) — x(^M(a) ~ ^M(a)) 



ia(s + s*)/2 


—i + a(s — s*)/2 


cosfca; 

(s — s*)/2 — ia 
—i{s + s*)/2 

+ sin kx 

—i{s + s*)/2 
— (s — s*)/2 — ia 


< 

i + a(s — s*)/2 


—ia{s + s*)/2 

> 


R M( t 3) - + ^M(/3)) 



2 + ip(t-t*) 


P(t + t*) 

cos kx 

t + t* 

2 p-i(t-t*) 

+ sin kx 

-i(t-t*)- 2 P 
-(t + t*) 


P (t + t*) 


2 -ip(t-t*) 





ip(t +1*) 


—2 i + Pit — t*) 

cos kx 

(t-t*)- 2 iP 
—i(t +1*) 

+ sin kx 

-i(t + t*) 
-{t-t*)-2ip 


2 i + P(t — t*) 


-iP(t +1*) 


Now let us detail the squaring method relied in the Majorana basis. Starting with 


(*7 M^a + M) 


di + M d t - d 3 0 -d 2 

—dt — c?3 —d\ + M d 2 0 

0 d 2 d±+ M —d t — d 3 

—d 2 0 d t — 83 —d\ + M 


and taking <f> = e ^ e lkiX e lk ^y e lk 3 Z ; we derive an explicit form for a matrix of four solutions 


{vlv} = e ~ ikx 


iki + M —ie — ik 3 0 

ie — ik 3 —iki + M ik 2 

0 ik 2 iki + M 

—ik 2 0 —ie — ik 3 


—ik 2 

0 

ie — ik 3 
—ik\ + M 


( 5 . 5 ) 
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This matrix can be decomposed into real and imaginary parts, The real part reads 

{Rj\ = 


M cos kx + ki sin kx —e sin kx — k 3 sin kx 0 -k -2 sin kx 

e sin kx — k 3 sin kx M cos kx — fci sin kx ^2 sin kx 0 

0 &2 sin kx M cos kx + fci sin kx e sin kx — k^, sin kx 

—fc2sinfcx 0 —esinfcx —/C3 sinfcx M cos kx — k\ sin kx 


Because the determinant of this matrix does not vanish 

det{i?j} = ( M 2 cos 2 kx — (fc 2 + fc 2 — e 2 + k 2 ) sin 2 kx) 2 = M 4 ; 

its columns represent linearly independent solutions of the Majorana equation. 

In similar manner we consider the matrix {Ij}: 


[-0] — * 


fci cos kx — M sin kx —e cos kx — k^, cos kx 
e cos kx — k$ cos kx —k\ cos kx — M sin kx 
0 k 2 cos kx 

—fc’2Cosfcx 0 


0 —fc2Cosfcx 

k 2 cos kx 0 

fci cos kx — M sin kx e cos kx — k$ cos kx 
—e cos kx — cos kx —k\ cos kx — M sin kx 


its determinant is 


det{/j} = ((fc 2 + fc| — e 2 + k 2 ) cos 2 kx — M 2 sin 2 kx) 2 = M 4 , 


therefore its columns dive linearly independent four solutions of the Majorana equation (with opposite 
charge parity). 

Let us write down explicit form of four real solutions: 


Ri = 


i?3 = 


and four imaginary solutions 


h=i 


h = i 


M cos kx + fci sin kx 


—e sin kx — sin kx 

e sin kx — fc’3 sin kx 

, i?2 — 

M cos kx — fci sin kx 

0 

k 2 sin kx 

—fc2 sin kx 


0 

0 


—fc’2 sin kx 

fc2 sin kx 

, i ?4 = 

0 

M cos kx + fci sin kx 

e sin kx — fc3 sin kx 

—e sin kx — fc3 sin kx 


M cos kx — fci sin kx 

ations 



fci cos kx — M sin kx 


—e cos kx — fc3 cos kx 

e cos kx — fc’3 cos kx 

0 

II 

<SJ. 

—fci cos kx — M sin kx 
k 2 cos kx 

— fc2 cos kx 


0 

0 


— fc’2 cos kx 

fc2 cos kx 

, I4 = i 

0 

fci cos kx — M sin kx 

e cos kx — k% cos kx 

—e cos kx — ks cos kx 


—fci cos kx — M sin kx 


( 5 . 6 ) 


( 5 . 7 ) 


Now, let us specify squared solutions in Majorana basis, when starting with the real scalar function 
(kx = et — kx)): 


$ = cos kx, (ilM^a + M) cos kx = {$j} 


fci sin kx + M cos kx 
e sin kx — k 3 sin kx 
0 

—fc’2 sin kx 


—e sin kx — fc’3 sin kx 
—fci sin kx + M cos kx 
k 2 sin kx 
0 


0 — fc2sinfcx 

fc2 sin kx 0 

fci sin kx + M cos kx e sin kx — fc’3 sin kx 
—e sin kx — k% sin kx —fci sin kx + M cos kx 


( 5 . 8 ) 
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note identity {4>j} = {Rj} ~ see (15.61) . 
Similarly, one derives 


fci cos kx — M sin kx 
e cos kx — k 3 cos kx 
0 

— k 2 cos kx 


<&' = —i sin kx, 


—e cos kx — k 3 cos kx 
—k\ cos kx — M sin kx 
k 2 cos kx 
0 


(ilMda + M) sin kx = {$'■} 


0 

k 2 cos kx 

fci cos kx — M sin kx 
—e cos kx — k 3 cos kx 


— k 2 cos kx 
0 

e cos kx — k 3 cos kx ’ 

—fci cos kx — M sin fcx 

(5.9) 


again note identity {$'} = {Ij}- 

It is readily verified, that these two 4-dimensional spaces, Rj and Ij, cannot be connected by any 
linear transformation. To this end, one should check the relationships 


hi) — '(j)kR(k) 


hi) 1 — S(j)kR(k)i ■ 


(5.10) 


Indeed, from (15.101) we infer that the matrix S depends on coordinates 


I = S R 


S = I R~ l 


i 

M 2 


x 


Mk\ — F sin2fcx 
— (e — k 3 ) 2 sin2fcx 


0 


— k 2 sin2fcx 


— (e + k 3 ) 2 sin2fcx 0 

— (Mfci + F sin 2kx) —k 2 sm2kx 

— k 2 sin2kx Mk\ — \{M 2 + k 2 ) sin2fcx 

0 — (e +/C3) 2 sin2fcx 


— k 2 sin2fcx 

0 

— (e — k 3 ) 2 sin2fcx ’ 
— (Mfci + F sin2fcx) 


where for brevity we use notation \{M 2 + fc 2 ) = F. It means that these two sets, {Ij} and {-Rj}, 
determine substantially different spaces: from physical standpoint they refer respectively to Majorana 
fermions with different properties with respect to the charge conjugation. 

In the charged Dirac case, at fixed momentum there exist sets of two states different in helicity, 
whereas in Majorana real case we face sets of four solutions only. 


6 Vanishing the current J z on the boundaries of the domain 
between two planes 


Using Dirac plane waves with fixed polarization states (for brevity the general factor e Kt e zklX e lk2V is 
omitted): 



1 


1 


k\+ik,2 


ki+ik2 

= * a ,k = e ikz 

k+p 

a 

, *2 = 4 h P) ,k = e ikz 

k—p 

P 


k 1 +ik 2 

k+p 


0 k 1 +ik 2 
^ k—p 


and similar ones t with the change k => — k: 



1 


1 


ki+ik2 


fcl -\-ik 2 

*3 = *a,-k = e~ lkz 

k—p 

a 

, ^4 = ®0) I -* = e ikz 

k+p 

P 


a ki +ik 2 

K—p 


0 k\+ik2 
r k+p 


( 6 . 1 ) 


( 6 . 2 ) 


let us make a linear combination $ = + A 2 ty 2 + ^4.34'3 + ^4^/4 . Four components of this wave 


function $ are 

$1 = e ~ ut e iax [ e lkz (A 1 + A 2 ) + e~ ikz (A 3 + A+\ , 
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$ 2 =e~ iet e iax (k 1 +ik 2 ) 


Jkz / A 1 


e lkz {- 


— e 


—ikz 


^3 


'fc+p k—p k—p k+p 

$ 3 = e - iet e“ x [e ifez (A ia + A 2 /3) + e~ ikz (A 3 a + A 4 p)] , 


$4 = e~ iet e iax (k\ + ik 2 ) 


Akz ( A\Oi . 2/3 \ — ikz/ ^-3^ . ^4/3 


e lkz (- 


-) — e~ lkz ( 


k—p k+p 


'k+p k — p 
Note the structure of the current (in spinor basis) is 

J z = $ + 7 ° 7 3 <i> = - <&3$ 3 ) - ( < f’ 2 < i ) 2 ^ $4$4) . 

The current on the boundaries (z = —a, +a) vanishes if the following requirements are fulfilled: 


(6.3) 


+ 1 = e ip $i 


$4 = e lCT $ 2 


[e~ lka (A ia + Ail 3) + e lka {A 3 a + A+ 3)] = e ip [e~ zka (A 1 + A 2 ) + e ika (A 3 + M)\ , 


„-ika, A l a , ^2/3 


A 3 a A 4 /3 




'k+p k—p 


'k—p k+p' 


ika 


A\ 


k+p k—p 


^2 ) _ e ifea ( ^3 


a 4 


k—p k+p 


z = +a , $ 3 = e ip ’$ 1 , $4 = e iv <&2 


[e lfca (Aia + A 2 p) + e~ ika (A 3 a + A 4 p)\ = e* [e ika (A 1 + A 2 ) + e~ ika (A 3 + A 4 )\ 


Aia A 2 /3 


A^q. A\ 


^ a (P^z + - e ~ ika (T^- + T~T ~~) 


'k+p k—p 


'k—p k + p' 


= e 


e ika { 


A, 


A 2 


k+p k—p 


) - e~ lka { 


A, 


^3 


k—p k+p 


Thus, we arrive at the homogeneous linear system with respect to complex variables A±, A 2 , A 3 , A 4 . 
It is convenient to introduce notation K = e 2 “ fe , then the above system reads as 

A ± (a - e ip ) + A 2 {p - e lp ) + A 3 (a - e ip )K + A 4 {P - e ip )K = 0 , 

ArK^a - e ip ) + A 2 K{p - e ip ) + A 3 (a - e ip ) + A 4 {P - e ip ) = 0 , 

Ai(a - e lcr )(k - p) + A 2 (p - e ia )(k +p) - A 3 K(a - e ia )(k +p) - A 4 K{P - e ia )(k — p) = 0 , 

A 4 K{a - e iv ){k - p) + A 2 K{p - e iv ){k+p) -A 3 {a- e iv ){k + p) - A 4 {P - e iv ){k - p) = 0 . 


(6.4) 


Let us write down an explicit form of the main matrix for linear system ( 16 . 41 ) 

(a - e ip ) (/3 - e ip ) (a - e ip )K (p - e ip )K 

(ot-e ip )K ( p-e ip )K (a - e ip ) {p - e ip ) 

(a — e ia )(k + p) {P — e ia )(k — p) —(a — e llT ){k — p)K —{p — e ia ){k+p)K 
{a — e lu ){k + p)K (/? — e lv ){k — p)K —(a — e ll ')(k—p) —(P — e lu )(k + p) 


(6.5) 


Eq. det S = 0 reduces to a 4-th order polynomial with respect to K, from physical point of view we 
are interested in roots which are complex number K = e 2lka of the unit length - see the detailed analysis 
of its solutions in [5j. 
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7 Vanishing the current on the boundaries, basis produced by 
the squaring method 

Now let us derive explicit form of vanishing the current J z on the boundaries of the domain, when using 
the set of solutions obtained within the squaring method. We start with (the general factor e — * et e lfc l a: e* fc2^, 
is omitted) 


= 


M sin kz 


0 

0 

, ^2 = 

M sin kz 

(e sin kz + ik cos kz) 

— (fci — ik 2 ) sinfcz 

— (fci + ik 2 ) sinfcz 


(e sin kz — ik cos kz) 



(e sin kz — ik cos kz) 


(fci — ik 2 ) sin kz 

^3 = 

(fci + ik 2 ) sin kz 

M sin kz 

, ^4 = 

(e sin kz + ik cos kz) 

0 


0 


M sin kz 


four components of the linear combination $ = Ai'S’i + A 2 ^ 2 + A3^3 + are 

$1 = A\M sin kz + A2O + Aa(esin kz — ik cos kz ) + A 3 (k\ — ik 2) sin kz , 
$2 = AiO + A 2 Msinkz + A 3 (ki + ik 2 ) sinfcz + ^(esinfcz + ikeoskz) , 
$3 = Ai(e sin kz + ikeoskz ) — A 2 {k\ — ik 2 ) sinfc^ + A 3 M sinfcz + A4O , 
<f>4 = —Ai(fci + ik 2 ) sinfcz + ^(esinfcz — ikeoskz) + A3O + A 3 M sinfcz . 

Remembering on the current structure 

J z = $ + 7 ° 7 3 $ = ($*$1 - $3^3) - ($2^2 - ^4^4) ■ 

we impose restrictions: 

z = -a , $3 = e v $i, $4 = e ltT (!> 2 => 

Ai(e sin ak — ik cos ak) — A.2(fci — ik 2 ) sin ak + A 3 M sin ak 
= e lp [AiM sin ak + A3 (e sinafc + ik cos ak) + ^4(^1 — ik 2 ) sinafc] , 

—Ai(fci + ik 2 ) sin ak + ^(esinafc + ik cos ak) + A4M sinafc 
< 5 > 2 [A 2 M sin ak + + ik 2 ) sin ak + ^(esinafc — ikeosak )] ; 


= e irT 


(7.1) 


(7.2) 


z = +a , $ 3 =e ip $i, 


$ 4 = A" 




Ai(esinafc + ikeoska) — A 2 (k\ — ik 2 ) sinafc + A 3 M sinafc 
= e zp [AiM sin ak + A 3 (e sinafc — ikeosak) + A4(fcq — ik 2 ) sinafc] , 
—Ai(fci + ik 2 ) sinafc + A2(e sinafc — ikeosak) + A 3 M sinafc 
= e lu $ 2 \A 2 M sinafc + ^(fcq + ik 2 ) sinafc + A^esinafc + ifccosafc)] . 


With the use of elementary identities 


n iak 


n —iak 


cos ak = 


K + K~ 


sin ak = —i 


giak _ g—iak 


K — K~ 


(7.3) 


and the shortening notations 

e + fc = m, e — fc = n , fci + ik 2 = / , fci — ik 2 = g , mn = fg , g = f* , 
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we arrive at linear equations 


= y, 


e — w 


— Ai[K(m — xM) — (n — xM)\ + A 2 g{K — 1) — A 3 [K(M — xn) — (M — xm)} + A±f? p g(K — 1) = 0 , 
—A\[K{n — yM) — (m — yM)] + A 2 g{K — 1) — A 3 [K(M — ym) — (M — yn)\ + A 4 yg{K — 1) = 0 , 
Aif(K — 1) — A 2 [K{n — vM) — (m — vM)} + A 3 vf(K — 1) — Ai\K(M — vn ) — (M — vm)\ = 0, 
A x f{K - 1) - A 2 [K (to - e^M) - (n - e iv M)} + A 3 e iv f{K - 1) - A 4 [I<{M - e i!/ m) — (M — e^n)] = 0. 


This linear homogeneous system has solutions if its determinant equals to zero: 


S' = 


(7.4) 


—K(m — xM) + n — xM 
—K(n — yM) + (m — yM ) 

f{K - 1) 
f(K- 1) 


g(K - 1) 

9(K - 1 ) 

—K[n — vM) + m — vM 
—K(m — wM) + n — wM 


—K(M — xn) + M — xm 
—K(M — ym) + M — yn 
vf(K - 1) 
wf(K - 1) 


xg(K — 1) 

yg(K - 1) 

—K(M — vn) + M — vm 
—K(M — wm) + M — wn 


In the paper ]6j it was shown that any particular choice in 4-dimensional space of states for the Dirac particle 
does nod influence substantially the result of solving the equations det S = 0 and det S' = 0: the whole set of the 
roots is the same. 


8 2-component Weyl neutrino 

The possible way to specify neutrino wave equation is to consider the Dirac equation in spinor basis 


i/j(x) = 


r](x) 


.£(*) = 


e 1 

£ 2 


g(x) = 


Vi 

V2 


,7 


0 <j a 

a a 0 


where a a = (/, +a k ) , cr a = (I, —cr k ) 


1 0 

1 

0 1 

2 

0 -i 

3 

1 0 

0 1 

, <7 = 

1 0 

» <7 = 

+i 0 

» <7 = 

0 -1 


Un this basis, we have two equations 

i<j a d a £(x) = Mg(x), ia aJ d a g{x) = M£(x), 
setting here m = 0 we produce equation for neutrino and anti-neutrino: 

(neutrino) i a a d a g(x) = 0 , (anti-neutrino) i a a d a £(*) = 0 . 
The neutrino equation differently reads 

(■ id t - idja j )ri = 0, E y = — id t g , 
where Pjtj 3 = E stands for helicity operator. For the plane wave 


— iet ikx 

7] = e e 


VI 

V2 


from (HOJ it follows that helicity of the neutrino with fixed momentum is automatically negative: 

E ?7 = — eg , e = + \jkj + k% + fc| . 

For those states one finds restriction on the components gi, 772 in (18.51 : 


e + k 3 

ki — ik 2 

Vi 

k x + ik 2 

e- k 3 

V2 


n _ . 2 ,2 . , 2 . , 2 fei + ik 2 

= 0 => e = k x + k 2 + k 3 , g 2 = -:— gi 

e — k 3 


( 8 . 1 ) 


( 8 . 2 ) 

(8.3) 

(8.4) 

(8.5) 


( 8 . 6 ) 
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Let the neutrino plane waves be written in the form (for simplicity, 771 = 1; and /C 3 = k; we will need also 
solutions with opposite k) 


— iet ixk 1 iyko izk 

77 = e e L e y z e 


1 

/ —iet ixk-t iyko —izk 

fj — g g e e 

1 

k\+ik2 


fci+ifc 2 

e — k 


e+/e 


(8.7) 


Summing equations for 3 and r] + : 

V + (dt- dj <J 3 )v = 0 , V + (dt - dja J )rj = 0 
one derives expression for the neutrino current 

d a J a = 0 , J a = (v + v, — r] + a j r i'J . 

Its components are 

J* = mvi + V 2 V 2 , J z = -vlvi + V 2 V 2 , 
J 1 = -vlv2 - vim , J 2 = ivlvi - ivlvi ■ 

In particular, the component j z is given by 

k 


( 8 . 8 ) 


3 = ~ 1 - 


fc? + kl 

(e - fc ) 2 


e — k 


(8.9) 


Now, let us examine the vanishing of the current J z on the boundaries of the domain between two planes. 
The structure of the current indicates the way to reach this: 


J z = -7]* ri 1 + rj 2 V 2 


V2 = e^Vi ■ 


( 8 . 10 ) 


To satisfy this requirement, let us introduce a special linear combinations of the plane waves with opposite 
momentums (general factor is omitted) 


H = Ar/ + Br[ = 


Ae izk + Be~ izk 

_k 1 ±ik2 Ae izk _ k 1 ±jk2 Be -izk 


for which the above current restrictions read 

z = —a : 

-Ae 


H 2 (—a) = e ip Hi(—a) 

Be iak = e ip [Ae~ iak + Be™*] ; 


kl + \ _ — iak ^1 + j-, _iak _ ip r a —iak , o ^iak 


t — k t + k 

z = +a : 7 / 2 ( 0 ) = e a H\(a) , 

k 1 + . iak k\ + ik2 j-. —iak ia r a iak . o —iak-\ 

t—k t + k 

Let us introduce parameter e 2lafc = K, then we get the linear homogeneous system with respect to A an B\ 

A{e ip + kl + f 2 ) + BK(e ip + kl + f 2 ) = 0 , 

e — k " ' u 


e + k 


AK(e ia + kl + f 2 ) + B(e ia + fcl + f 2 ) = 0 . 


e — k 


With the shortening notations 

ki + ik 2 
e — k 


= f 


ki + ik2 
t + k 


t + k 


= 9, /<?* = !, /*<?=!; 


the system reads 


A(e ip + f)+BK ( e ip + g ) =0, 
AK (e ia + f)+B [e ia + 3 ) = 0 . 


( 8 . 11 ) 
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Nontrivial solution exists if K 2 satisfies the equation below 


K 2 (e^ + f)(e^+g) 
(e ia + f)(e i P + g) ' 

Let us simplify the notation: e lp = x , e ICT = y, then 

K 2 = (x + f)(y + g) . 
(x + g)(y + f) ’ 


( 8 . 12 ) 


(8.13) 


It is ready proved that the right hand part of (18.1311 indeed is a complex number of the unit length. To this end, 
it suffices to multiply it by the complex conjugate - so we obtain the needed identity 

(x + f)(y + g) (1 /x + f*)(l/y + g*) _ (x + f)(y + g) ( g + x) xf (y + f) yg _ 1 

(x + g)(y + f) (l/x + g*)(i/y + /*) (x + g)(y + f) ( f + x)xg (y + g)yf~ 

then eq. (18.131) reads as the quantization rule for k in the form 


g 4 iak _ (e lp + f)(e la + g) 

~ (e ip + g){e™ + f) ' 


(8.14) 


This equation permits different solutions. Mostly of them can be found only numerically. For the case e lp = e ICT , 
we have an analytically solvable equation, e 4lafc = 1. 


9 On the representation of J z = 0 in arbitrary basis of the Dirac 
matrices 


In spinor basis, vanishing of the current is achieved if 

4/3 = e' p ’F 1 , 4/4 = e ICT 4' 2 ; 


(9.1) 


however these relations are not covariant: they are valid only in this particular basis. This deficiency can be 
resolved. Let us present m in the matrix form $ = G$: 


4-r 


o 

T 

O 

O 


4-r 

4*2 


0 0 0 e~ irT 


4> 2 

4< 3 


O 

o 

o 

y 

% 


4'3 

4>4 


0 e iCT 0 0 


4' 4 


The matrix G can be decomposed in terms of Dirac matrices in spinor representation - such form will be auto¬ 
matically covariant (valid in all bases) 


0 a a 

5 

o 

1 

5 a 

0 — (j a 

a a 0 

5 7 = 

0 I 

,7 7 = 

<r“ 0 


Decomposition for G (19.21) should be searched in the form 

G = (n 0 7° + n 3 7 3 ) + 7 5 (m 0 7° + m 3 ■y 3 ) ; 


which is equivalent to 


G = 


0 

n 0 + m 0 + (n 3 + m 3 )<T 3 


no — mo — (n 3 — m 3 )(j 3 
0 


(9.3) 


and leads to the system 

[n 0 - mo) - (n 3 - m 3 ) = e~ zp , 

(n 0 — mo) + (n 3 - m 3 ) = e~ z<T , 

(no + m 0 ) + (n 3 + m 3 ) = e +zp , 

(no + mo) — (no + m 3 ) = e + * CT . (9.4) 
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Its solution is 


1 . . 1 ,. . . . , 
n 0 = — (cos p + cos a) , ri3 = - (* sm p — * sin a) , 

1 ,. . ... 1 , 

mo = -(jsmp + jsma) , m3 = - (cos p — cos a) . 


Correspondingly, the matrix G is presented as follows 

G = i(cos P + cos < 7 ) 7 ° + i(i sinp — isincr) 7 3 
+7 J [i (i sin p + i sin < 7 ) 7 ° + i (cos p — cos er) 7 3 ] ; 


This can be translated to other form 
G = e 


_ , ip 1 + 7° 7 U +7 J ^ 1-7° 7 U ~7 J 


+ e 


, +ia 1 + 7 5 7°-7 3 , —itr 1 -7 5 7 0 + 7 3 

+ e --— +e —2-— • 

Condition for vanishing the current in covariant form is given by the relation 

J z = 0 G(p, < 7 ) S'= $ . 

Formula (ESP becomes simpler at e ICT = e lp : 


G = e 

Similarly, (19.51) is simplified if e ICT = e ~ lp : 

G = e 


_ a +i P 1 + 7 5 „,o 0 -i P 1 - 7 5 . 


7 +e“ 


_ „+»p 7° + 7 3 t -i„ 7°~7 3 


(9.5) 


(9.6) 


(9.7) 


(9.8) 


10. Conclusion 

The accent in the paper was given to consideration of solutions of the Dirac equation which have vanishing the 
third projection of the conserved current J z on the boundaries of the domain between two parallel planes. Such 
solutions are reachable when considering 4-dimensional basis of four solutions - plane waves with opposite signs 
of the third projection of momentum +£3 and — £ 3 . It is shown that the known method to solve the Dirac 
equation trough squaring method, if based on the scalar function 4? = e~ x,Lt e xklx e xk2V sm(kz + a), leads to 4- 
dimensional basis of Dirac solutions. It is shown that so constructed basis is equivalent to the 4-space of Dirac 
plane wave states; corresponding linear transformation is found. Application of the squaring method in Majorana 
representation is investigated as well. 

General conditions for vanishing third projection of the current J z at the boundaries of the domain between 
two parallel planes are formulated: on the 4-dimensional base of plane spinor waves and on the base of the squared 
basis. In both cases, these conditions reduce to linear homogeneous algebraic systems, equation det S = 0 for 
which turns to be a 4-th order polynomial, the roots of which are e 2lka , where a is a half-distance between the 
planes. Solutions of this 4-th order polynomial will be considered in a separate paper. 

The case of Weyl neutrino field is investigated as well; conditions for vanishing the neutrino current J z on the 
two boundaries reduce to 2 -nd order polynomial: its solutions are found; all of them are complex numbers with 
the unit length. 
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